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Abstract If massive sterile neutrinos exist, their de- 
cays into photons and/or electron-positron pairs may 
give rise to observable consequences. We consider the 
possibility that MeV sterile neutrino decays lead to the 
diffuse positron annihilation line in the Milky Way cen- 
ter, and we thus obtain bounds on the sterile neutrino 



decay rate > 10 



-28 



from relevant astrophysi- 



cal/cosmological data. Also, we expect a soft gamma 
flux of 1.2 X 10"'* - 9.7 X 10"'' ph cm'^ s'^ from the 
Milky Way center which shows up as a small MeV bump 
in the background photon spectrum. Furthermore, we 
estimate the flux of active neutrinos produced by ster- 
ile neutrino decays to be 0.02 
through the earth. 
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1 Introduction 

Understanding the nature of dark matter remains a 
fundamental problem in astrophysics and cosmology. 
Since the discovery of neutrinos' non-z ero rest mass 
(jFukuda et al.lll998t iBilenkv et al.lllQQSl) . the possibil- 
ity that neutrinos contribute to cosmological dark mat- 
ter has become a hot topic again. In particular, the ster- 
ile neutrinos belong to a class of candidate dark matter 
particles with no standard model interaction. Although 
the recent MiniBooNE data challenges the LSND result 
that suggests the existence o f eV scale sterile neutrinos 
( Aguilar-Arevalo et al.ll2007 ) , more massive sterile neu- 
trinos (eg. keV, MeV) may still exist. The fact that ac- 
tive neutrinos have mass implies that right-handed neu- 
trinos should exist which may indeed be massive ster- 
ile neutrinos. The existence of the sterile neutrinos has 

M. H. Chan and M. -C. Chu 

Department of Physics and Institute of Theoretical Physics, 
The Chinese University of Hong Kong, 
Shatin, New Territories, Hong Kong, China 



been invok ed to explain many phenomena such as miss- 
ing mass (|Dodelson and Widrowl Il994t IShi and Fuller 
1999) and the high temperature of the hot gas in Milky 
Way and clusters (jChan and Chull2007l[20oi ). There- 
fore, it is worthwhile to discuss observational conse- 
quences if massive sterile neutrinos exist, which may 
decay into light neutrinos, positron-electron pairs and 
photons. In this article, we consider the possibility that 
sterile neutrino decays give rise to the 511 keV lines in 
Milky Way and thus obtain bounds on the mass rus and 
total decay rate F of the sterile neutrinos using relevant 
observational data. 



2 511 keV photon flux 



The bright 511 keV annihilation line from Milky Way 



has been observed for a few decades ( Leventhal et al 



19781 : Knodlseder et al. 2005 ). and its origin has been 
much debated. Recent values of the 511 keV pho- 
ton flux from the bulge and disk are (1.05 ± 0.06) x 
10-3 ph cm-' ^ s"^ and (0.7 ± 4) x l O'^ ph cm'^ s'^ 
respectively ( Knodlseder et al. 20051) . Assuming a 
positronium fraction of fp = 0.93, one can trans- 
late these intensities to annihilation rates of (1.5 ± 
0.1) X 10'^^ s~^ arid (0.3 ±0.2) x 10^^ g-i respectively 



(jKnodlseder et al.l 120051 ). The annihilation rate in the 
bulge is several times larger than that in the disk. 
The source of positrons in the disk can be explained 
by the decay of ^^Al. Using a disk model, the pho- 
ton flux is calculated to be 5 x 10"'* ph cm~^ s~^, 
which can account for 6 — 100% of the disk emission 
20051 ). However, the origin of the 



( Knodlseder et al. 



bulge source is still an open question in astrophysics. 
Recent observation by INTEGRAL/IBIS indicates that 
the upper limit of photon flux for reso lved single point 



sourc es is 1.6 x 10 ph cm s (|De Cesare et al 



20051) . which means that the 511 keV annihilation line 



2 



comes from mainly diffuse sources ratlier than point 
sources. 

There have been many models trying to explain 
the 511 keV annihilation line of bulge emission. Po- 
tential sources include neutron stars or blackholes 



2.2 Positron channel 

The positrons produced will be slowed down due to 
ioni zation losses. The power loss is approximately given 
by (|Longaiiill98ll ) 



Lingenfelter an d Ramaty 19837 , supernova remnants 



Dermer and Murphv.2001.). Wolf-Rayet stars (Ramatv et a^m ^ _2 x 10^^" f - ) (ln7 + 6.6) eV/; 



19791). pulsar wind dChi et al.lll996HWang et al.ll2006l) . 
and Gamma Ray Bursts ([Lingenfelter and Hueteijll984l) .| 

None of these can provide a satisf actory explana t ion as 



(5) 



Boehm et all 



they are mainly point sources, 
proposed the anni hilation of dark matter as a d iffuse 
source. Recent l y. IPicciotto and Pospelovl (j2005l ) and 



Khalil and Setol ( 2008 ) suggested that heavy sterile 
neutrinos with nis > 1 MeV can be a diffuse source 
of the 511 keV photon flux. In the following, we extend 
the idea from these two papers and obtain bounds on 
the decay rate of sterile neutrinos from observational 
data of the 511 keV annihilation line. 

2.1 Decay of Sterile Neutrinos 

A sterile neutrino i>s can decay into an electron-positron 
pair, photons and lighter neutrinos ly through different 
chan nels. The major c hannel is I'g with decay 

rate ( Barger et al.lll995l ) 



where 7 is the Lorentz factor of the positron and n is 
the average number density of electrons in the galactic 
bulge. The stopping distance d for 1 MeV positrons in 
this process is about 10^^ cm. Also we should consider 
the magnetic field in the Milky Way. The Larmor radius 
of positrons with energy is given by 



eB 



10 



13 



104 MeV / V 10-5 G 



B 



(6) 



The magnetic field strength in the Milky Way is about 
B — IQ-^ G. Therefore, for a 1 MeV positron, the 
Larmor radius is about 10^ cm. The stopping distance 
for the simple random walk of a positron, the distance 
that a positron is co nfined, is about \/rd ^ 1 pc or less 
( Boehm et al.l 120041 ) , which is much shorter than the 
mean free path of the annihilation: 



(1) 

where Gp and 9 are the Fermi constant and mixing 
angle of sterile neutrino with active neutrinos respec- 
tivel y. The radiatiye cha nnel is — > + 7 with decay 
rate ([Barger et al.lll995l ) 



/g± = ^ 30 kpc, 

n(Ta 

"vBiere 



meC^(7 -|- 1) 



7^-1-47-1-1 



r 



1 



ln(7 + A/72 - 1) 



(7) 



7 + 3 
(8) 



10247r 



7r4 ^ VlkeV/ 

(2) 



is the cross s ection of electron-positron annihilation 
jHe itler 1954). Therefore, the positrons will become 
-iilii-relativistic before annihilation. However, the rate 
fol a positron to annihilate with an electron in the dif- 
fuse region of Milky Way is 



where a is the fine structure constant. An electron- 



positron pair is produced through Vs ^ v + e"' 
with decay rate (jPicciotto and Pospelovll2005l ) 



P ^ nOaVe ~ 10 



-18 



G\ 

3847r3 



1 



\v\ 



1 



(9) 



0.1 cm-3 (jMuno et al.l 120041 ) and 



where n 

10^ cm s~^ is the mean sp eed of electrons in Milky Way 



where \V\ < 1 is a parameter, 
decay rate is 



(3) 

Therefore, the total 



r 



3iy 



fMv\ 



V4|yp 



(4) 



( Marconi and Hu"ntll2003 ). In order to produce lO^'^ s~^ 
annihilations, there must exist a large positron cloud 
with lO^'^ positrons in the Milky Way, and the initial 
production rate should be much greater than the anni- 
hilation rate. 

Suppose a sterile neutrino halo is formed and the 
positron production rate is much higher than the an- 
nihilation rate during the galaxy formation due to the 
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small n in the protogalaxy. The positrons will accu- 
mulate in the protogalaxy. The rate of change in the 
positron number density is given by 



n„+ 



n„+n(TaVe 



(10) 



where ns{t) is the number density of sterile neutrinos 
in the Milky Way at time t. Since we have ns{t) — 
nsoe~^^, where Uso is the initial number density of ster- 
ile neutrinos, the solution of Eq. (10) is 



'>''dt' + C 



(11) 



where C is a constant. After a long time, assuming an 
equilibrium is established at present time to, we have 
ns(tQ)Te = nf,+ naaVe- The total annihilation rate in 
the bulge is given by 



A 



bulge 



(12) 



where Rb is the bulge radius, which is assumed to be 
2.40 — 3.71 kpc in the model used bv iKnodlseder et al. 
(j2005r ). Similarly, the annihilation rate in the disk is 
given by 



A, 



disk ' 



4:7Trhns{to)Tedr, 



(13) 



where Rd and h are the radius and half of the thickness 
of the disk respe c tively . In the disk models used in 
Knodlseder et al.l (|2005[ ). the maximum Rd is 15 kpc, 
and the scale heights of young and old disk models are 
70 pc and 200 pc respectively. Here we assume that the 
sterile neutrino profile follows the dark matter profile, 
which can be mo delled by the isothermal n(s) = n nr~^ 
or NFW profile (|Navarro. Frenk and Whit(j|l996l) : 



r/a{l + rjaf 



(14) 



where a and nj, are parameters in the NFW profile. In 
the isothermal profile, the ratio of the annihilation rate 
in the bulge to that in the disk is 



A 



bulge 



Rb 



6-13. 



Adisk h\xY{RDlRB) 

In the NFW profile, the ratio is 



(15) 



A 



bulge 



idisk 



In 



Rb 



1 



Rsja 
RB/a + l 



RbIo. 



- 1 
(16) 



Since over 60% of disk emission can be explained by the 
decay of ^^Al, the lower bound of the ratio of the dif- 
fuse emission should be about 7. Therefore, the isother- 
mal profile agrees better with the observed ratio. Fur- 
thermore, since ^buigo = (1-5 ± 0.1) x lO"*"^ s^^, from 
Eq. (12), we have n^Y ^ - 10^^ m^^ s'^ The up- 
per limit of central mass density in isothermal model 
constrains mgriQ < 5.5 x 10^^ kg m^^. For nis > 1 
MeV, we get no < 3 x 10^^ m'^. Therefore, F > 3 x 
10~^^ s~^. Similarly, for NFW model, we have npFe ~ 
4 X IQ-^^ m-^ s-i and m.nf, < 3.2 x 10^22 ^^-3^ 
For rus > 1 MeV, we get Ug < 2 x 10* m~'^ and 



Fe > 2x10" 



IfF, 



10-^« s-^ andm^ ~ 1 MeV, 



we can get sin 29 ~ 10 , which is co nsistent with the 



diffus e X-ray background constraint (jBovarskv et al 
20091) . 



2.3 Radiative channel 

There exists another radiative decay channel which 
gives a photon flux with energy ma/ 2: 



sFds. 



(17) 



line of sight 



The b ranching ratio is given by (jPicciotto and Pospelov 
20051) 



0.031 

4|y|2-Hi' 



(18) 



and therefore the photon flux should be = (1.2 x 
10"" - 9.7 X 10"^) ph cm-2 g-^. Basically, aU the 
emitted photons are monochromatic with energy E = 
ms/2. However, some of the photons will scatter with 
intersteller medium before reaching us. The probability 
of the scattering is 



'line of sight 



(19) 



where dc is the Compton cross section. For n w 1 cm^'^, 
the total Pg from the disk and bulge is w 1 x 10^^. Due 
to the scattering, the energy distribution of the photons 
reaching us is broadened slightly. Fig. 1 shows the con- 
tribution of the emitted photon flux (we assumed E — 1 
MeV) together with the diffuse background photon flux 
dF/dE = 2 62(^/ 0.1 MeV)-2-^5 MeV"! cm-^ 



(Kinzer et al 



19971 ). The emitted photons contribute 
2 —15 ,9cLpf the background flux at around 1 MeV 
— which s Hows?8 am4ll'MeV bump' in the spectrum. The 
ump is a clcisical problem in observational as- 
tronomy which is long conjectured to be a real feature in 



the spectrum (jKinzer et al.l 119971 ). However, the MeV 



bump is now commonly believed to be an artifact of 
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incomplete background rejection (jKinzer et al.lll997l ). 
Nevertheless, at least in this model, part of the MeV 
bump is a real feature in the diffuse background photon 
spectrum. 

2.4 Active neutrino channel 

The lighter active neutrinos are produced in the main 
decay channel. The total active neutrino flux due to 
sterile neutrino decays passing through the earth is 
0.02 — 0.1 cm~^ s~^. The total number of active neu- 
trinos passing through IceCube - the larges t neutrino 
detec tor in the world - is about 10^ s~^ (jWiebusch 
20091) . Although this flux is theoretically detectable, 
the energy of the decayed neutri nos is too small to be 
detected in current experiments (jLunardinil 120061 ) . 

Active neutrinos may interact with neutrons or pro- 
tons in a pulsar to produce electrons or positrons. The 



3 Discussion and Summary 

The fact that active neutrinos have finite masses implies 
that right-handed neutrinos should exist which may in- 
deed be massive sterile neutrinos. The existence of the 
sterile neutrinos has been invoked to explain many as- 
trophysical ph enomena such as the cooling flow prob- 
lem in clusters ( Chan and Chj2007l) . In this article, we 
consider the possibility that the decays of MeV sterile 
neutrinos act as a source of the 511 keV flux line. The 
decaying sterile neutrinos provide a diffuse source of 
positrons which can account for the required flux. The 
large bulge to disk ratio of 511 keV luminosity can also 
be accounted for if the decaying sterile neutrinos follow 
the isothermal distribution. From the observed 511 keV 
photon flux in the Milky Way, we obtain the allowed 
ranges of the sterile neutrino decay rate F > 10~^^ s~^. 



cross section of such interactions is (T^ - 10-^\EjlOMeY)^}k¥'^^ we do not have representative tight bounds 



where E^, is the energy of the neutrinos. For example, 
in a typical pulsar, the average number density is about 
10^® cm~^, and the mean free path for a 1 MeV neutrino 
in the pulsar is 10^ cm. As the crust thickness is also 
of order 10^ cm, almost every neutrino passing through 
a pulsar will interact with the neutrons and protons to 
produce electrons and positrons. As a result, a huge 
amount of electrons and positrons is produced and af- 
fected by the strong magnetic field B in the pulsars to 
emit synchrotron radiation. The synchrotron frequency 
of the electrons is given by 



2'iTnipC 



= 2.8 X IQ^^-il 



B 



1012 G 



Hz, 



(20) 



where 7e is the Lorentz factor of the electrons. There- 
fore, the frequency of the synchrotron radiation lies 
within the x-ray band. The power emitted by one elec- 
tron is 



2„4o2 2 



B 



1012 G 



erg s 



on decay rate of sterile neutrinos, the results are still 
compatible with cosmologi cal bounds and cluste r cool- 
ing fiow (F < 10-1^ s-i) (|Ghan and Chull2007t) . The 
radiative decay channel produces soft gamma rays, with 
an expected flux of 1.2 x 10"''- 9.7 x 10"'' ph cm"2 s"i, 
which show up as a small MeV bump in the back- 
ground photon spectrum. The total active neutrino 
flux due to sterile neutrino decays is estimated to be 
0.02 — 0.1 cm~2 g-i in the vicinity of the earth. These 
active neutrinos interact with neutrons and protons in 
pulsars to produce x-ray photons which may be de- 
tectable in the future. 
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(21) 



For a pulsar with radius 10 km nearby, the total num- 
ber of active neutrinos due to sterile neutrino decays 
passing through the pulsar is IO12 s~i. Assuming 
B = IO12 G and 7e = 4, the total power emitted is 
^ 10^^ erg s"^, which is much less than the upper limit 
of the non-therm al x-ray luminosity in a typical pul- 



sar (10^° ergs 1) ( Zavlin and Pavlovl2004 ) . Therefore 



only a very small peak near the synchrotron frequency 
may appear in the x-ray spectrum. Assuming that the 
10^ or so pulsars in Milky Way all see similar neutrino 
flux, the resulting synchrotron radiation ca n contribute 



abou t 102'' erg s 1 to the background x-ray (jLvne et al, 
19981) . 
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Fig. 1 The spectrum of the background photons including the MeV photons coming from sterile neutrino decays. 



